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Abstract 

For a system composed of two particles Bell's theorem asserts that averages of physical quantities 
determined from local variables must conform to a family of inequalities. In this work we show 
that a classical model containing a local probabilistic interaction in the measurement process can 
lead to a violation of the Bell inequalities. We first introduce two-particle phase-space distributions 
in classical mechanics constructed to be the analogs of quantum mechanical angular momentum 
eigenstates. These distributions are then employed in four schemes characterized by different types 
of detectors measuring the angular momenta. When the model includes an interaction between the 
detector and the measured particle leading to ensemble dependencies, the relevant Bell inequalities 
are violated if total angular momentum is required to be conserved. The violation is explained 
by identifying assumptions made in the derivation of Bell's theorem that are not fulfilled by the 
model. These assumptions will be argued to be too restrictive to see in the violation of the Bell 
inequalities a faithful signature of nonlocality. 

PACS numbers: 03.65.Ud,03.65.Ta,45.20.dc 
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I. INTRODUCTION 

Bell's theorem was originally introduced [l, 2| to examine quantitatively the consequences 
of postulating hidden variable distributions on the incompleteness of quantum mechanics 
put forward by Einstein, Podolsky and Rosen [sl (EPR). In particular, the hidden variables 
were supposed to locally and causally complete quantum mechanics by making sense of the 
'reality' of physical quantities described by non-commuting operators relative to two spatially 
separated particles in an entangled state. Bell showed that a correlation function obtained 
from averages over the hidden variables of these physical quantities must satisfy certain 
inequalities (the Bell inequalities), and that these inequalities are violated by quantum 
mechanical averages. Given that experiments have confirmed with increasing precision the 
correctness of the quantum formalism, it is generally stated that the violation of the Bell 
inequalities contradicts locality. The strong version of such statements asserts that quantum 
mechanics itself is non-local [4]. This vocable is quite popular (in particular among non- 
specialists as well as in quantum information papers) but there is a general agreement 

nriQ 

among most specialists that this strong assertion is unsubstantiated [5|, la, lZ|] . Instead, the 
received view is the weak version following which Bell's theorem asserts the incompatibility 
of local hidden variables with quantum mechanics. Nevertheless it can objected, in principle 
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or through abstract models 
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whether the assumptions made in order to derive 



Bell's theorem are necessary in order to enforce locality, or whether they only rule out a 
certain manner of ascribing local variables to the measurement outcomes. 

In this work we will show that statistical distributions in classical mechanics can vio- 
late Bell-type inequalities. Moreover the statistical distributions we will employ are not 
exotic objects but the classical analogues of the quantum-mechanical coupled angular mo- 
menta eigenstates, so that our model is essentially the classical version of the paradigmatic 
2-particles singlet state. The violation of the inequalities can of course be achieved only 
provided the model falls outside the assumptions necessary in order to prove Bell's theo- 
rem. This role will be played by a probabilistic interaction that is assumed to take place 
between the measured particle and the detector, combined with the requirement that the 
total angular momentum be conserved. Although this interaction is local, it nevertheless 
spoils the derivation of Bell's theorem, because it introduces an ensemble dependency of the 
outcomes: the resulting averages then involve correlations given by conditional probabilities 
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between ensembles rather than between the individual phase-space positions. As a conse- 
quence, the different expectation values employed in Bell's inequalities cannot be derived 
jointly, as required in the derivation of the theorem. 

The paper is organized as follows. Wc will start by introducing the classical phase-space 
distributions (Sec. 11), first for a single particle, then for 2 particles with total zero angular 
momentum. We will explain why these distributions are the classical analogues of the 
quantum angular momentum eigenstates. In Sec. Ill we will investigate three different 
examples of Bell-type models. Each of the examples will be characterized by the same phase- 
space distribution but by differing detection schemes. In the first case, the projections of the 
angular momentum of each of the particles along arbitrary axes are directly measured by the 
detectors, leading to a straightforward application of Bell's theorem (which will be briefly 
derived). In the second example the detectors yield discrete outcomes, depending on the 
values of the angular momenta; this example, which also abides by Bell's theorem, will allow 
us to introduce conditional probabilities to account for the correlated angular momenta. 
The third example will illustrate the same situation with stochastic variables (the angular 
momenta specify probabilities of obtaining an outcome). In Sec. IV, we will introduce 
an example falling outside the class of Bell-type models. This example will also involve 
discrete measurement outcomes, but the presence of an interaction leading to ensemble 
dependencies will be introduced. We will see that ensemble dependencies lead to non- 
commutative measurements for a single particle, and to the violation of the Bell inequalities 
for initially correlated two-particle systems. In Sec. V we will discuss these results, insisting 
on the role played by the existence of joint distributions and on the relationship between 
locality and conservation laws. A short summary and our conclusion are given in Sec. VI. 



II. CLASSICAL DISTRIBUTIONS ANALOGUES OF ANGULAR MOMENTA 
EIGENSTATES 

A. One-particle angular momentum distributions 

A quantum mechanical angular momentum eigenstate Ijoui) is characterized by a well- 
defined value a/ Jo (jo + 1) of the modulus of the angular momentum J and of its projection 
(of value m) along a given axis z. In configuration space the spherical harmonic | {6, jom) f 
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gives the probability distribution corresponding to a fixed value of J and Jz as 6 and (the 
polar and azimuthal angles) span the unit sphere. The classical statistical distributions can 
be considered either in phase-space, defined hy VL = {0, 4>,pe,P4,} where pg and p^j, are the 
conjugate canonical momenta, or in configuration space. Let us assume the modulus J of 
the angular momentum is fixed. Let Pz{^) be the distribution in phase-space given by 

Pz,ie,<p,pg,p^) = NS{JM - JzMJ'm - J',). (1) 

defines a distribution in which every particle has an angular momentum with the same 
magnitude, namely Jq, and the same projection on the z axis Jz^, without any additional 
constraint. Hence pz^ can be considered as a classical analog of the quantum mechanical 
density matrix Ijo"^) Oo"^!- Eq. ([T]) can be integrated over the conjugate momenta to yield 
the configuration space distribution 



sin{e)Jji-Jl/sin\e) 
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(2) 



where we have used the defining relations J^{^) = p^ and J^{Q) = pi ^p^^j sin^ Q to perform 
the integration. Further integrating over Q and and requiring the phase-space integration 
of p to be unity allows to set the normalization constant N = Jo/27r^. 

p{6, (p) gives the statistical distribution of the particles in configuration space. Its stan- 
dard graphical representation (parameterization on the unit sphere) is shown in Fig. 1(a) 
along with the quantum mechanical orbital momentum eigenstate (a spherical harmonic 
taken for the same values of j and m) in Fig. 1(b). The similarity of both figures is a state- 



ment of the quantum-classical correspondence in the semiclassical regime, since \J p{6, 0) is 
approximately the amplitude of the configuration space quantum mechanical eigenstate for 
high quantum numbers. Note that rather than working with the particle distributions in 
configuration space, it will also be convenient to visualize the distribution of the angular 
momentum in physical space corresponding to a given particle distribution (see Fig. 1(c)); 
6 and (p will then denote the position of J on the the angular momentum sphere. 

Let us take a second axis a making an angle 6a relative to the z axis (in this paper we will 
take all the axes to lie in the zy plane). We can define a distribution by fixing the projection 
Ja of the angular momentum on a to be constant, pa^ = S{Ja — Jao)S{J — Jq) ■ In configuration 
space, this distribution may be shown to be obtained by rotating the distribution of Eq. ([2]) 
by the angle 6a towards the a axis. We will be interested below in determining the average 




FIG. 1: Normalized angular distribution for a single particle in configuration space, (a) Classi- 
cal distribution p{6,(p) of Eq. ([2]). (b) Corresponding quantum distribution (spherical harmonic 
\YjM{0,(j))\^ with J/r] = 40, r/ = h, and M/J = 5/8 as in (a)), (c) Distribution of the angular 
momentum on the sphere for a distribution of the type p{9,(j)), invariant around the z axis with a 
fixed value of Jz- 

projection Ja on the a axis for a distribution of the type ([2]) corresponding to a well defined 
value of Jz. Using Ja = Jz cos 9a + Jy sin 9a, Jz = P(j> and 

{Jy)j^ = J J sin 9 sin (f) sin 9 aS{p(j, — J zQ)dQ = (3) 

by rotational invariance, we obtain 

{Ja)j^^ = j P4> COS 9a5{p^ - Jzo)dVt = J z^ COS 9a. (4) 

Note that a given J can belong jointly to several distributions pa,, and pf,g {a and b being 
different directions). But if we require that any distribution must correspond to a well- 
defined value of the angular momentum projection along a given axis, then distributions 
such as pa^ and become mutually exclusive. The ring-like distribution of the angular 
momentum on the angular momentum sphere represented in Fig. 1(c) (corresponding to the 
configuration space distribution shown in Fig. 1(a)) can be generalized to cover the entire 
hemisphere centered on the z axis (see Fig. 3(a)). Then properties such as Jia and Jib being 
of the same sign on such hemispheres become mutually exclusive properties. 
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B. Two-particle angular momentum distributions 



The situation we will consider below, by analogy with the well-known EPR-Bohm pairs 
in quantum mechanics, is that of the fragmentation of an initial particle with a total angular 
momentum Jj- = into 2 particles carrying angular momenta Ji and J2. Conservation of 
the total angular momentum imposes Ji = J2 = J and 



Eq. ([5]) implies a correlation, imposed initially at the source, between the angular momenta 
of the 2 particles and of their projections along any axis a: the knowledge of the value of 
Jia allows to infer the value of J2a, J2a = —Jia- Without further constraints (or additional 
knowledge), the classical distribution in the 2-particle phase space is given by 



where is again a normalization constant. The corresponding distributions of the angular 
momenta in physical space - easier to visualize than p - is uniform on the sphere, with Ji 
and J2 pointing in opposite directions (see Fig. 2(a)), reflecting the isotropic character of 
the fragmentation as well as the correlation ([5]). 

III. BELL-TYPE MODELS 
A. Setting 

The Bell inequalities are obtained by computing average values of measurement outcomes 
performed independently on each of the 2 particles. Three examples are studied below, all 
involving the initial fragmentation of a particle with zero angular momentum (Sec. II. B). In 
the first example, we assume that the measurements give directly the value of the projection 
of the angular momentum of each particle along an arbitrarily chosen axis. In the sec- 
ond example we introduce detectors having a threshold, resulting in discrete measurement 
outcomes depending solely on the position of the particles' angular momenta. The third 
example is a repetition of the second but with stochastic variables. Bell's theorem, which 
is derived in Sec. III.B, is verified in all these cases. To alight the notation, we will choose 
units such that J = 1. 



Ji + J2 = 0. 



(5) 



p{ni, ^2) = N5iJi + J2)5( Ji - J'), 



(6) 
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B. Bell's theorem 



1. Example 1: direct measurement of the classical angular momenta 

Two particles with initial total angular momentum Jt- = flow apart. Let a and h be 
two axes in the zy plane. The projection of particle I's angular momentum along the a axis, 
Jia and that of particle 2 along 6, are measured. The average of the joint measurement 
outcomes on the 2 particles is directly given by the values of Jia and J2b and the probability 
distribution given by Eq. ([6]). All these quantities depend on the phase-space position of 
the particles, i.e. on the position of the angular momenta on the sphere (see Fig. 2(a)). The 
average is computed from 

{JlaJ2b) = j JU^l)J2b{^2)p{^l,^2)dnidQ2. (7) 

Given the rotational symmetry, z is chosen along a, hence Jia = and 

■hh = P,p2 cos {6b - 6a) + { J sin 62 sin 02 sin {6b - 6a)} ■ (8) 

One first integrates over 02 (the term between {..} vanishes) then over p^^ (yielding p^^ = 
—p^i because of the correlation 5 (Jia + J2a))- The last non-trivial integration is over p^-^. 



{JiaJ2b) = j dp^, - pI^ cos {6b - 6a) 2tiN j dVL 



(9) 



where dVt represents the variables remaining after the integration of the delta functions. 
Since p is normalized, we have 

j dp^^27iN j dn = l. (10) 

Integrating Eq. (fTOj) over p^-^ allows to obtain the value between the [..] in Eq. IQ thereby 
avoiding the explicit calculation of the normalization constant. The result for the expectation 
is 

E{a, b) = {JiaJ2b) = -I COS {6b - 6a) . (11) 
2. Derivation of the Bell inequality 

The correlation function C{a, b, a', b') involved in Bell's inequality is given by 

C(a, 6, a', b') = {\E{a, b) - E{a, b') \ + \E{a\ b) + E{a\ b')\) /V^^ (12) 



where a' and h' are arbitrary axes in the xy plane and Knax is the maximal absolute value 
that can be obtained in a measurement outcome. Let us denote by Via(f2i), V2b{S^2) etc. the 
detected values along the relevant axes, with the 2-particle average being 

E{a,h) = j VU^^)V2h{^2)p{^i.^2)dn^dn2. (13) 

The Bell inequality 

C(a,6,a',6') < 2 (14) 



puts a bound on the value of the correlation function. It is obtained 12| by forming the 
difference 

E{a,h) - E{a,b') = J Vi„(fii) [^26(^^2) - V2b'in2)] p{n,n2)dn,dn2 (is) 

where Via has been factored. Likewise, 

E{a', b) + E{a', h') = J Via' [V^26 + V^26'] pdnidn2. (16) 
We now use IV2/3I < Knax {13 = b, b') to derive 

I V2b - V2y \ + \V2b + V2V I < 2\4.ax. (17) 

Take the absolute values and use |Viq,| < V^ax (« = a, a') in each of the Eqs. ffTSl) and 
(|T6l) to obtain two inequalities. Adding these inequalities and using (fTTll leads to the Bell 
inequality 

\E{a,b) - E{a,b')\ + \E{a' ,b) + E{a' ,b')\ < 21^1,. (18) 

In the present example, Knax = 1, and C{a,b, a' ,b') is bounded by 2^/2/3, so that the Bell 
inequahty f|T4l) is verified. 

As a corollary, note that the factorization made in Eqs. (fT5|) -( IT6|) is equivalent |l^ to 
the existence of joint distributions of the form 

:Faba'b' = J Via{ni)V2b{n2)Via'{ni)V2b'{n2)p{nun2)dnidn2. (19) 



Indeed, Bell's inequality can be proved |14] by adding and substracting J^aba'b' from Eq. ([T5|) 
and then factorizing ViaV2b and ViaV2b' respectively. The term Taba'v is the average obtained 
when 4 measurements are made - 2 outcomes are obtained for each particle (particle's 1 
V property is measured along the axes a and a! whereas particle 2 is measured along the 
axes b and b'\ The factorization, or equivalently the existence of Taba'b'-, is an important 
assumption in the derivation of the inequalities. 
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3. Derivation in the stochastic case and joint distributions 



In the stochastic case, a given phase-space position (f^i, Q2) does not determine a unique 
valued outcome (Via(fii), V2b{^2)) as above (corresponding to what is usually termed "deter- 
ministic case") but determines instead well-defined probabilities piVia, V26, fli, ^^2) of obtain- 
ing (Via, V2b)- The counterpart of the factorization made in Eq. (ITS!) hes in the factorization 
of the probabilities, 

p{Vla, V2b, ^2) = p{Via, n,)p{V2b, ^2), (20) 

where p{Via, fli) is the single particle elementary probability such that 

P{Via) = j p{v,,,ni)p{ni)dQ,. (21) 

The expectation value (fT3i) is then replaced by 

E{a,b) = j VU^i)V2b{^2)p{^u^2)dnidn2 (22) 

with 

VU^^) = Y,yiaP{Vla,^l) (23) 
V2b{^2) = Y.'y^bP{V2bM2). (24) 

The derivation leading to Eq. f|T8l) proceeds as above by replacing the value of the outcomes 
by their respective averages Via and V2b- The factorization (l20l) allows to obtain a joint 
probability for an arbitrary number of events from the elementary probabilities p(V, fi); the 
counterpart to Eq. (IT^ is 



aba'b' 



piy^a. ^l)piy2b. ^2)piyia'.^l)piy2b'. ^2)P{S^U ^2)dnidn2. (25) 



Note that the existence of a joint probability Paba'b' (that appears here as a consequence 
of the factorization (!20|) ) leads immediately to the inequality (ITSl) irrespective of any other 
assumption concerning the dependence of the outcomes or probabilities on supplementary 
variables (here the phase-space positions, the 'hidden- variables' in quantum mechanics). 
Indeed, using expressions of the type 

E{a,h)= J2 ^1-^26 J2 ^-ba'b' (26) 
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for the average values, we have 

\E{a, h) - E{a, 6') I < Paba'b' I V^la {V2b - V2b')\ (27) 

and an analog inequality for \E{a', h) + E{a\ h')\. Adding both inequalities yields 

\E{a, h) - E{a, h') \ + \E{a\ h) + E{a\ h') \ < 

J2 P-^-'"' (1^1- (^26 - V2b')\ + \Via' {V2b + V2b')\) < 2V^,,, (28) 

since the expression between (...) is bounded by 2V^^^ and the joint probability sums to 1. 
C. Discrete outcomes 

In this second example, we take over the setup of the first example except for the mea- 
surement outcomes: we now assume that a given detector placed on an axis can only give 
2 values, depending on the sign of the angular momentum's projection. The outcomes are 
given by 

f i if Ji„ > f i if J26 > 

D.^m = {\ D2bm = {\ (29) 

y -\ if Ji„ < y -\ if J2b < 

and depend only on the positions Ji and J2 of the angular momentum (hence on the phase- 
space position of the measured particles). The average value 

{D^aD2b) = j DU^,)D2t{n2)p{^i,^2)dnidQ2 (30) 

takes the form 

1/2 

{DiaD2b) = / Pd^idn2 (31) 

where 'D{k, k') represents the domain of integration on which the joint conditions sign( J2a) = 
— sign(Jia) = — sign(/c) and sign(J2b) = sign(/c') hold (see Fig. 2(b)). The integral gives the 
probability 

Pkk' = P{D^a = kr\D2b = k') = P{D,a = k)P{D2b = k'\D,a = k) (32) 

where P{D2b = k'\Dia = k) is the probability of obtaining D2b = k' conditioned on the 
knowledge that Dia = k. The conditional probability appears because of the initial correla- 
tion ([5]) - the positions of the angular momenta are not independent. The conditional proba- 
bility can more easily be determined on the angular momentum sphere by computing the area 
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where sign( J2fe) = sign(fc') relative to the area of the hemisphere where sign( J2a) = — sign(fc) 
(of area 27r). This area is given by the intersection of the two relevant hemispheres (see Fig. 
2(b)), i.e. a spherical lune whose area can be put under the form 2'Kk{k — k') +4kk'{6b — da)- 
Since p is uniform on the sphere, we have P{Dia = k) = 1/2 from where 

2kk' 

Pkk' = k{k - k') + — \e,-ea\, (33) 

TT 

and the average {DiaD2b) becomes 

E(a,6) = -i + ^^. (34) 

The maximal detected value here is Knax = 1/2- The correlation function is computed from 
Eq. ( |T2|) and it may be verified that C(a, 6, a', h') is bounded by 2: Bell's inequality is 
again verified. 



D. Discrete outcomes: a stochastic model 



We now elaborate on the preceding example to give a model in line with the stochastic 
version of Bell-type variables. A given position of the angular momentum of a particle in 
phase-space does not specify the outcome S, as in Eq. fl29l) . but the probabilities p{Sia = 
k, Qi) of obtaining the outcome k. For definiteness we will replace Eqs. (!29l) by 

1 f |if > 1 f |if J26 > 

2 1 r ^ n 2 11..^ 





1 Uif Jia > 1 iif J2b > 

p{S,a = --,n^) = {l p{S2b = --, ^2) = <^ ; . (36) 



4 if J2b < 



2 [fifJi„<0 2 

The expectation value involves first averaging, for each phase space position, over the two 
possible outcomes, before averaging over the distribution p of the angular momenta: 

{SiaS2b) = j su^i)S2bmp{nu^2)dnidQ2 (37) 

with (cf. Eqs. (I22])-(I24D) 

Siam = J2kp{Sia = k,ni) (38) 

k 

S2b{^2) = ^'Pi^^b = k', Q2). (39) 
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Taking into account the correlation at the source [Eqs. ([5])- ([6])], we proceed as in the 
preceding example, except that now each probability Pkk' contains several contributions 
with a weight given by p{Sia = k, Qi)p{S2b = k', Q2) that depends, through Eqs. (jSSD-dSSD, 
on the domains ©(±1/2, ±1/2) over which sign(Jia) = =1=1 and sign(J26) = =tl. For example 
for k, k' = ^, we have 

Pi 1 = — / pdQidQ2-\ / pdQidQ2-\ / pdQidQ2-\ / pdVtidVt2] 



16^(1,1) 16^(i,_i) 16^(_i,i) 16jc(_i,_i) 

(40) 

now each integral represents a probability P(sign(Jia) = =f1 n sign(J2fe) = ±1). Comparing 
with Eqs. fl5^ - fl52]) . we see that in the stochastic case, the probabilities Pkk' depend as 
in the preceding example on the areas on the angular momentum sphere occupied by the 
individual positions of each angular momentum compatible with the outcomes (although in 
the stochastic case there are many more such areas, each contributing with a given weight). 
Overall, Eq. ([37D yields 

{S,aS2b)= ^^'^''' = \{^^^^-^- 

k,k'=-l/2 ^ ^ 

C{a, a', b, b') is readily computed and is again, in line with Bell's theorem, bounded by 2. 



IV. A DETECTION MODEL VIOLATING THE INEQUALITIES 

The fourth example has similarities and differences with the models studied in Sees. 
III.C and III.D. A given detector on an axis measures the angular momentum's projection 
but only delivers the outcomes ±1/2. However, the outcomes depend on a probabilistic 
random interaction between the detected particle and the detector. This interaction has a 
specific property (it vanishes on average) that results in the introduction of an ensemble 
dependency. We will see that this feature combined with the conservation of the angular 
momentum between ensembles prevents the factorization that was seen above to be necessary 
in order to derive Bell's theorem. 
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FIG. 2: (a) Uniform distribution of Ji and J2 on the unit sphere; the angular momenta are 
correlated via the conservation law ([5]) and must thus point in opposite directions. In the first 
example (Sec. III.B), the detectors measure Jia and the correlated J2b as the angular momenta 
span the sphere, (b) Example 2 (Sec. III.C): Distribution of J2 , when Dia = — 1/2 was obtained. 
A measurement of D26 will yield ±1/2 depending on the position of J2: if J 2 lies within the 
light shaded region (intersection of the two positive hemispheres centered on a and on b, denoted 
P(— 1/2, 1/2) in the text), D26 = 1/2 will be obtained, —1/2 when J2 belongs to the dark-shaded 
region (2?(-l/2, -1/2)). 

A. Particle-detector interaction for a single particle 

1. Basic properties 

Let pi(r2i) be the phase-space distribution for the single particle 1 and Ria = ±1/2 
denote the outcome obtained by placing a detector on the a axis. Let P{Ria = k, pi) be the 
probability of obtaining the reading k on the detector if the statistical distribution of particle 
1 (or equivalently, the distribution of Ji) is known to be pi. We will impose the following 
constraint on the interaction: the average {Jia) over phase-space of the measured value is 
the one obtained by averaging over the measurement outcomes. This constraint takes the 
form 

1/2 

(^i«)p. = E kP{R,, = k,p,) = {J,,)^^, (41) 

fc=-l/2 

meaning that whereas individual outcomes depend on the interaction, on average the net 
effect of this interaction is zero. The models leading to Eq. (HT!) are not unique - any model 
verifying Eq. pTl) and obeying Ylk^^^^a ~ k,Pi) = 1 will do. Depending on the specific 
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model, Eq. (HT!) will not be verified for an arbitrarily chosen distribution pi; only a class of 
distributions can be consistent within a given model. In the present model, we will assume 
as in the previous examples that pi can only be a uniform distribution occupying one (or 
both) of the two hemispheres of the angular momentum sphere. 

Let us examine for such distributions the consequences of Eq. (HTj) . Assume that pi 
corresponds to a uniform distribution of J i on the positive hemisphere centered on the a 
axis, to be denoted pia+- If a measurement is made along the h axis, a direct computation 
of (^ib)pi.+ gives 

A; 

If one measures Ria the average (H2l) becomes +1/2, i.e. the only positive detected outcome. 
Therefore, since the probabilities are positive, we must have 

P(i?i, = l/2,pi«+) = l (43) 
P(i?i, = -l/2,pi,+) = 0. (44) 

Conversely if the distribution is pia- (uniform on the lower hemisphere) we obtain the 
opposite probabilities 

P(i?i, = l/2,pi,_) = (45) 
P(Pia = -l/2,pi„_) = l (46) 

and 

{Rih)p,,_ = J2 ^P^J^^ = ^' P^-) = -\ COS - ■ (47) 

k 

Note that Eq. (H2l) along with the normalization of the probabilities uniquely determines 
the value of the probabilities 

n/n I 1 \ C0s(6'b -6a)±l . 

P{Rib = = ^ ^ (48) 

as well as the equality between the relative expectation value corresponding to positive 
(resp. negative) outcomes Pi;, and the average of the angular momentum projection over 
the regions where Jn, is positive (resp. negative), i.e. 

± ^P(Pi, = pi„+) = {H{±Ju)J^,)^^^^ , (49) 
14 



H denoting the unit-step function. 

The main property of this particle-detector interaction based model is that the detected 
result does not depend on a phase-space point (or on a given individual position of the 
particle's angular momentum on the sphere), be it through a deterministic value ascription 
or through well-defined probabilities. Indeed, if this were the case, then Eqs. fl43l)- (HB|l 
would imply that 

R,ai^i) = l/2^J,a>0 and = -1/2 ^ Ji„ < 0, (50) 

as in the example studied in Sec. III.D. But then assume Rib is measured and the ensemble 
is known to be pia+ (uniform distribution on the positive hemisphere centered on the a axis). 
On the angular momentum sphere pia+ can be seen as being composed of the intersections 
with pib+ and pu-, Pia+ = (Pia+ ^ Pib+) U (pia+ npi5_). The respective integration domains 
are V{—^, i) and V{—-^, — |) (we use the notation introduced in Sec. III.C; see Fig. 2(b))). 
Hence 

4. = / ^ib(^i)pia+(^i)c^^i = ^ / pib+{Qi)dQi - ]- I pib-{Qi)dni, 

(51) 

yielding (1 — 2{9b — 9a) /t^)/2 in contradiction with the constraint ( HTl) defining the model, 
{M,,^^ = cos {9b - 9a) 12. 

We see therefore that the value ascription given by Eq. fISUl) does not fit with the main 
property of the model ^. The reason is that Eq. fHTl) introduces an ensemble dependency on 
the model: the probabilities do not depend on the phase-space position but on the ensemble, 
as if the particle's angular momentum effectively occupied an entire hemisphere (physically, 
this may happen for example if the particle follows a stochastic motion with its angular 
momentum constrained to remain in the ensemble, the timescale of the measurement being 
significantly larger than the timescale of the stochastic motion). Eq. fl50l) should thus be 
replaced by 

Riai^i) = ±1/2 ^ Ji„ ^ for every Jia e pi. (52) 



Since value ascriptions given by V{Vl) or p{V, Vl) are characteristic of Bell models, it could be said that 
even for a single particle, Bell-type models are inconsistent with the present model. On the other hand, 
it could be argued that the ensemble pi should be taken as the 'hidden variable', given that this is the 
variable ascribing values to the outcomes and probabilities, even though pi may not qualify following 
Bell's terminology as being a 'beable' (see Sect. V). 
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FIG. 3: (a) The ensemble pia+ for the single particle model described in Sec. IV. B. Any Ji in 
this ensemble has a positive projection Jia, measuring Ria gives the outcome +1/2 with certainty, 
without changing the ensemble, since pia+ is symmetric relative to the a axis and {Jia) p^^^ = 1/2. 
(b) In the same situation Rib is measured. Now the symmetry axis of the ensemble pia+ does not 
coincide with the b axis. Hence measuring Ri^ can yield either +1/2 or —1/2 with probabilities 
proportional to {H{zizJib)Jib)p^^^- The ensemble is modified during the measurement, undergoing 
a rotation toward the positive or negative b axis as indicated by the arrows, (c) The two-particle 
distribution after Ria was measured and the outcome is known to have been Ria = —1/2, in which 
case particle 1 is described by the ensemble pia-- Conservation of the angular momentum then 
requires that particle 2 be described by P2a+, so that if i?2a were measured, the outcome +1/2 
would be obtained with certainty [Eq. (|58p ]. If instead i?2f) is measured, we have a single particle 
problem for particle 2 identical to the one portrayed in Fig. 2(b). 

2. Further considerations 

Although this has no effect on the computations, it will be convenient, in order to provide 
a physical interpretation, to detail the consequences arising from the model. Eq. fl52|) 
associates an outcome Ria along an axis a with Jia being of the same sign for every member 
of the hemispheric ensemble pia± (see Fig. 3(a)). Since Ria = ±1/2 = {Jia) p^^^, we can 
envisage that the random interaction occurring during a measurement effectively changes 
the distribution of the angular momentum: for instance if initially the distribution is on a 
given hemisphere, say pia+, Eq. fl52|) is realized and Ria = 1/2 is obtained with certainty, 
reflecting {Jia) p^^^ ■ If Rib is measured, the final distribution is pih+ (resp. pib-) if the 
outcome k = 1/2 (resp. —1/2) is obtained (see Fig. 3(b)). The outcome thus appears as the 
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average value of the angular momentum projection in the post-measurement distribution 
and Eq. fH2|) becomes 

A:=±l 

Note that this implies that consecutive measurements involving projections along different 
axes do not commute: the condition fl52l) cannot be realized jointly along two different 
directions, like the classical analogues of the angular momenta eigenstates presented in Sec. 
II ^. If the initial distribution is pia+ measuring Rn, then Ria' entails that Rn, is measured 
over pia+ but Ria' over pii,± depending on the outcome Rn,. In the reverse order, Ria' is 
measured first, the average being given by {Jic) p^^^ and Rih then involves the values of Jih 
over one of the distributions pia'±- Eq- fl53l) also allows to compute the change in the angular 
momentum projection due the measurement, 

A ( Ji,) = ( - ( = -2kP{R,, = k, pi,^.), k = (54) 

Consider now the uniform distribution on the entire sphere pis. It can first be envisaged as 
the angular momentum occupying the upper or lower hemispheres along a definite direction 
(say a) so that 

Pie = (pia+ + Pia-)/2. (55) 

Since distributions in classical mechanics obey the principle of linear superposition, pis can 
also be taken as a sum of the expressions given by Eq. fl55|) over different directions a. 
Alternatively the angle a in Eq. (155|) can be taken to vary in time (then the measurement 
does not involve a change in the distribution but rather a selection of the angular momenta 
such that Jia > or Jia < 0), or J be distributed on the entire spherical surface (then the 
measurement induces a change in the distribution pis Pia±)- Only in these latter cases is 
the distribution spherically symmetric; all these possibilities lead to the same probabilities 
and average values, yielding P{Ria = ±1/2, pi^,) = 1/2 for any axis a as well as a vanishing 
average fHTj) as required. 

^ We have chosen distributions on hemispheres rather than the ring distributions of Sec. II for continuity 
with the examples investigated in Sec. Ill; the model studied here would also hold if ring like distributions 
were employed. 
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B. Two-particle expectation 

1. Distribution and conservation of the angular momentum 

Before computing the two-particle averages and correlation functions, we explicitate the 
initial distribution and the conservation of the angular momentum for the model. We have 
seen that the defining property Eq. fHTl) implied that value ascription depended on distri- 
butions (taken to be ensembles on given hemispheres) and not on individual phase-space 
positions. The two-particle distribution given above by Eq. ([6]) is (i) spherically symmetric 
and (ii) anti-correlates the individual positions of the angular- momenta 3 2 = — Ji, so that 
we have = —Jia for projections along arbitrary axes a on the sphere. We require the ex- 
tension of these two properties so that they hold over the initial distribution, to be denoted 
by ps- We must thus have for any of the two particles i and axis direction a 

{R^a),^ = {J^a),^ = 0, (56) 

SO that both outcomes Ria = ±1/2 can be obtained with equal probability. The correlation 
between the outcomes for the two particles is obtained by applying Eq. (HOl) to ps_ giving 

{HiJ2a)J2a),, = {Hi-J^a)Jla),,- (57) 

By Eq. (!53|) we have Ria = {Jia) p^^^ so that by way of Eq. (!57j) the anti-correlation 
J2a = —Jia implies that the outcomes and the distributions for the particles along the same 
axis are anti-correlated, 

{J2a),,^^ ^ R2a = 'Ria ^ " (Jia),,,, • (58) 

Eqs. ( I56l) - (|58ll hold for any arbitrary axis a. The anti-correlation for the measurement 
outcomes, based on the conservation of the angular momentum over the ensembles, implies 
anti-correlations between these ensembles. Measuring Ria links the outcome to one of the 
two ensembles pia± depending on whether Ria = ±1/2. In turn, this also fixes p2 = P2a^- 
Note that contrarily to the correlation between individual phase-space positions (for which 
one has J2a = —Jia and J2b = —Jib for any axes a and 6), Eq. (!58|) cannot hold jointly 
along several directions. This is a consequence of Eq. ( l52l) not holding simultaneously along 
several axes. 
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There are different possibilities for choosing exphcit reahzations of ps : all these pos- 
sibilities lead to the same results and all hinge on the conservation of the total angular 
momentum along an arbitrary axis demanded by Eq. fl58p . For example ps can be taken 
as proportional to pib+P2b~+ Pib~P2b+- Eq- dSH]) is then ensured provided the change in the 
angular momentum after the first measurement is taken into account in the angular 
momentum balance for the other particle. Alternatively as in Eq. (155|) . b can be taken as 
varying in time, giving 

PS = ^ {plb(t)+P2b(t)- + Plbit)-P2b{t)+) ■ (59) 

As for the single particle distribution pis, measuring Ria then selects the individual positions 
of Jia such that Jia ^ 0, correlated to the individual positions J2a ^ 0. Another possibility 
for ps would be to take the distribution and consider Ria as inducing a change in the 
distribution ps Pia±- 



2. Computation of the correlation 



Since the measurement outcomes do not depend on the individual phase-space positions, 
the average E{a, h) = {RiaR2b) cannot be obtained as in the preceding example from the 
phase-space averages (l30|) -(l3T]). but from the probabilities of detecting a given outcome as 
a function of the distribution. E{a, b) is computed from the general formula, also employed 
in Sec. III.C 



{RlaR2b) 

where as in Eq. (!32|) Pkk' is given by 



1/2 

p^= ^ kk'Pkk' 

k,k'=-l/2 



(60) 



Pkk' = P{Ria = knR2b = k\ ps) = P{Rla = k)P{R2b = k'\Ria = k) 



(61) 



and the two particle expectation takes the form 



1/2 



{RlaR2b)p^= Yl kP{Rla = k) 

k=-l/2 



1/2 



k'P{R2b = k'\R,a = k) 

k'=-l/2 



(62) 



For any particle i and direction a, we have 



P{Ria = ±1/2, p^) = 1/2. 
19 



(63) 



The conditional probability P{R2b = k'\Ria = k) is, as in the example involving discrete 
outcomes studied above in Sec. IIII C\ the probability of obtaining R2b = k' if it known 
that Ria = k. But we have just seen that obtaining an outcome Ria = k is linked to the 
respective densities Pia[sign(fc)] and P2a[sign(-fc)]- The conditional probability is therefore given 
by 

P(i?26 = k'\Ria = k) = P{R2b = k',p2a[sign{-k)]), (64) 

which is a single particle probability of the type given by Eq. (148|) . Note that in order 
to compute the expectation value, we do not need to know the values of these individual 
probabilities, as the knowledge of the conditional expectation - the expression between 
brackets in Eq. (1^21) - is sufficient. The 2-particle conditional expectation is given by the 
single particle average (J2b)n whose expression was determined above (Eqs. (HTl) . 

P2a[sign( — fc)] 

(142|) and (147|) ). We can rewrite the average in the form 

1/2 

J2 k'P{R2b = k'\R^a = k) = -kcosieb-ea). (65) 

k'=-l/2 

We can now compute the expectation E{a,b) = {RiaR2b) from Eqs. flB^ and fl^ . The 
result is easily seen to be 

E{a,b) = -^cos{9,-9a). (66) 

In the present derivation, we have assumed that the knowledge of particle I's outcome was 
obtained first, hence the appearance of the conditional probability regarding the outcomes 
of particle 2. But obviously by Bayes' theorem the result is the same if we assume instead 
that R2b is known first, and the conditional probability then concerns the computation of 
the outcomes of particle 1. 

The correlation function C{a,b,a',b') is again given by Eq. (fT2|) with V^ax = 1/2- 
C{a,b,a',b') violates the Bell inequality f|T^ for a wide range of values, the maximal vi- 
olation being obtained for C(0, ^, ^,^) = 2-\/2. This correlation function, with E{a,b) 
given by Eq. (jHSD, is familiar from quantum mechanics - it is precisely the correlation ob- 
tained for the 2 particles with spin 1/2 in the singlet state. It was shown in this case that 
Eq. (166|) can be seen as a consequence of a particular correlation between vectors whose 



projection is conserved on average 



16|- 
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V. DISCUSSION 



A. Ensemble dependence 

We have seen in our fourth example (Sec. IV) that correlation functions obtained from 
2-particle distributions in classical mechanics can lead to a violation of the Bell inequalities, 
without nonlocality being explicitly involved (it may play an implicit role, see Sec. C below). 
The main difference between this model and the other examples we have given consists in 
the ensemble dependencies: probabilities, average values and conservation laws are relative 
to a collective property (a given distribution) and do not depend, as in the other cases, 
on the individual phase-space positions. Indeed, the constraint (HTi) cuts the link between 
a definite phase-space position of a particle and a given measurement outcome (be it in a 
probabilistic or deterministic way). 

In this respect, it is noteworthy to compare the interpretation of the conditional probabil- 
ities appearing in Examples 2 (Sec. III.C) and 4 (Sec. IV). In both cases P{V2b = k'\Via = k) 
is grounded on the correlation and represents the probability of obtaining = k' given 
the knowledge that Via = k. In both cases the distribution of J2 is modified once the 
outcome Via = k is known ^ (it changes from a uniform distribution on the sphere to a 
uniform distribution on the positive or negative hemisphere centered on a, depending on k). 
However in example 2 the probabilities depend on the individual phase-space positions of 
the particles: although it may be unknown in practice, Ji has in principle a definite position 
that determines Via = k, and to this position corresponds the definite position J2 = — Ji 
that will determine the outcome V2b', so the conditional probability is computed by finding 
the individual positions of J2 such that V2b = k' compatible with the positions of Ji imposed 
by Via = k (namely Jia > 0). In example 4 an outcome Via = k cannot be linked in principle 
to an individual position of Ji and thus we can only infer from the outcome the ensemble 
to which Ji must belong; then from the conservation law we know the distribution for J 2, 
which allows to compute P{V2b = k'\Via = k) from the probability P{R2b = fc', P2a[sign{-fc)])- 
Hence we can only correlate observable outcomes with ensembles, not with individual po- 

^ It seems it is necessary to stress that the change in the probabiUty distribution of particle 2 when the 
outcome of particle 1 is known is not a physical phenomenon involving action at a distance, but the result 
of the information brought by the knowledge of the first outcome, given the conservation law. This point, 
unrelated to Bell's theorem, is an elementary inference common in the calculus of probabilities. 
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sitions of the angular momenta. Assuming that a given phase-space position determines 
probabihties, as in the stochastic model of Sec. III.D, only brings in several combinations 
of possible outcomes allowed by the definite positions of Ji and J2 = — Ji on the angular 
momentum sphere, but still allows to correlate these individual positions with measurement 
outcomes. 

B. Joint distributions and non- commutative measurements 

We had remarked in Sec. III.B that the existence of a joint probability Paba'y is sufficient 
to ensure that a Bell-type inequality holds, irrespective of whether the assumption that 
measurement outcomes and probabilities depend on the individual phase space positions 
is made. But if that specific assumption is made, then one is lead to the existence of a 
joint probability by imposing the factorization fl20p . Along these lines, there are two ways 
of seeing why Bell's theorem does not apply to our fourth example. 

First, the ensemble dependence can formally be thought of as arising from elementary 
phase-space probability functions specific to a given ensemble, i.e. 



(compare with Eq. (pT!) ). By employing expressions such as Eq. (1671) in the expectation 
value as given by Eqs. fl60|) - fl6T]) . it can be seen directly that the ensemble dependence of 
the elementary probabilities spoils the factorization (120|) - for example instead of ^(^2); one 
has outcome dependent expressions such as ^(^^2; P2{Ria))- 

The second manner starts with the remark made above concerning the non-commutation 
of the R measurements introduced in our model; in classical mechanics, measurements usu- 
ally commute, but this is not the case if they arise from collective phenomena (encapsulated 
in the ensemble dependency). By requiring that the angular momentum be conserved be- 
tween ensembles (just as it is when individual positions are considered), the consequences of 
the non-commutation are carried over from one particle to the other. As seen in Sec. IIV A2I 
in the single particle case, the probabilities and outcomes when Rib is measured after a 
first measurement is made will be different depending on whether Ria or Ria' was measured 
first. Because Eq. ( !58|) links the outcomes with the ensembles, this is also the case in the 
two-particle problem when R2b (or -R26') is determined after Ria or Ria' were measured. Put 




(67) 
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differently, although Eq. (158|) holds for the a, a', b and b' axes, it cannot hold jointly for all 
the axes because the single particle ensembles pia± and pia'± are mutually exclusive, as well 
as p2b± and P2b'± (see Sec. II and Sec. IV. A. 2). Hence a joint probability Paha'w cannot be 
defined and the model is not constrained by the inequality (1251) . This is consistent with the 



equivalence 



13l . llTj l shown between the verification of Bell's inequality and the commuta- 



tion of the four observables entering Eq. (fT9l) . The ensemble-dependence introduced in our 
model appears as a tool in order to enforce, in a classical context, the non-commutation of 
the measurements along different axes made on the same particle. 



C. Conservation laws and locality 

Factorization, enforcing the existence of joint distributions, and as such a necessary as- 
sumption in the derivation of Bell's theorem, is usually argued to be intimately linked to 
locality. According to Bell [3], factorization is a consequence of local causality, given that 
space-like separated events can only have common causes in their backward light-cone: there- 
fore the probability of obtaining a certain result in an event regarding one of the particles 
cannot depend on what has been measured on the other. It is known however that factoriza- 
tion can be seen 



IS 



20| as a consequence of two separate conditions, outcome independence 
(the conditional probability of one event does not depend on the outcome obtained in the 
other event) and parameter independence (dependence on the measurement direction of the 
other event). Only the violation of outcome independence can result in a genuine violation of 
local causality (it would permit superluminal signalling), whereas the violation of parameter 



independence allows a 'peaceful coexistence' |20| between local causality and other types of 
correlations preventing the factorization. 

The present model - like many quantum mechanical entangled systems - respects parame- 
ter independence [Eq. fl^ ] but violates outcome independence [Eq. f lMl) ] (the dependencies 
here must be understood relative to the ensembles and not relative to the individual positions 
of the angular momenta). This outcome dependence of the conditional probabilities is due to 
the conservation of the angular momentum, as encapsulated by J2 = — Ji (anti-correlation 
between individual positions), Eq. (!57j) (correlation between ensembles occupying opposite 
hemispheres centered on the same arbitrary axis) and Eq. flSSl) (anti-correlation between 
the outcomes made on the same arbitrary axis). Parameter independence on the other hand 
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guarantees that the predictions relative to Ria do not depend on what measurement or 
whether a measurement is carried out on particle 2, and vice- versa [Eq. [63]) ]. It is clear nev- 
ertheless that the angular momentum conservation affects the distributions of both particles. 
For example if ps is given by Eq. flS^ . made up from rotating distributions, measuring Ria 
not only freezes the rotation of particle's 1 distribution, but also that of particle 2 (precisely 
because the angular momenta are correlated and need to be conserved). If ps is taken as 
a uniform distribution on the sphere for the individually anti-correlated angular momenta, 
measuring Ria changes the distribution pi^ pia± but also p2E — * P2a=F- Hence, it can 
be argued that the conservation of the angular momentum as implemented in our model 
actually results from an implicit implementation of nonlocality. There are several answers 
to this question, depending on the status one gives to conservations laws and ensemble dis- 
tributions, or on how nonlocality or causality are defined. The three following positions can 
be singled out: 

1. The changes in the distributions are real physical effects, but the conservation of the 
angular momentum results from a symmetry that is intrinsically linked to space-time. 
Indeed the correlation fISS]) arises by generalizing the angular momentum conservation 
for individual positions to ensembles accounting for non-commutative measurements. 
There is no need to invoke a specific mechanism for a conservation law - conservation 
laws and symmetry principles are just postulated. But if desired, a field can be can 
be ascribed the role of transporting the angular momentum; in this respect, it may 
be useful to make the analogy with Feynman's paradox in which mechanical angular 
momentum is transmitted between two charged particles through the electromagnetic 
field |21|. 



2. The changes in the distributions are real physical effects due to a nonlocal form of cau- 
sation. The requirement given by Eq. fl55]) is sufficient to imply nonlocality. Action at 
a distance effects are quite common in non-relativistic classical mechanics, although 
the modern view is to ascribe such effects (like gravity or several phenomena in electro- 
statics) to the action of fields. Here the nonlocal effect would consist in accounting for 
angular momentum conservation. This does not necessarily contradict the preceding 
position since it can be argued that symmetries can give rise to nonlocality, a position 
leading to a holistic vision of symmetries as holding beyond a space-time framework. 
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3. The changes in the distributions are not physical effects: one must distinguish the 
observed frequencies (which are measured) from the calculus of probabilities (whose 
role is to make logical inferences given a certain state of information I22I). Conditional 
probabilities do not therefore express causation and the factorization of the probabili- 
ties does not follow from the requirement of local causality. In Bell's term, the variables 
entering the probabilities are not beables, an argument that may be supported by the 
fact that individual angular momenta positions do not ascribe values and that the sta- 
tus of ensembles as beables is questionable. The ensembles and their correlations are 
theoretical constructs encapsulating the state of knowledge we have of the situation, 
including the constraints (like conservation laws). 



VI. CONCLUSION 



To summarize, we have first constructed classical distributions analogues of the quantum 
mechanical angular momentum eigenstates; these classical distributions are characterized 
by being mutually exclusive, leading, with appropriate assumptions to non-commutative 
measurements. We have then derived Bell's theorem in the deterministic and stochastic 
cases; both cases are characterized by the fact that an individual position of the angular 
momentum ascribes a value (with certainty or with a given probability) to a measurement 
of the projection along any axis. As a result, a joint probability distribution for an arbitrary 
number of events can be defined. Three different examples of Bell-type models were studied. 
A non-Bell-type model was introduced in Sec. IV: in this model, individual positions of the 
angular momenta are irrelevant to determine the measurement outcomes, that only depend 
on ensembles. As a result, single particle measurements do not commute, and a distribution 
for joint measurements along different axes cannot be defined. If it further assumed that 
the total angular momentum must be conserved, the Bell inequalities are violated. 

The present results do not disprove Bell's theorem - as we have seen, in these circum- 
stances the assumptions made in the derivation of the theorem are not fulfilled. We have 
argued that the violation of the inequalities in our classical model is due to the conservation 
of the total angular momentum in the context of non-commutative measurements; nonlo- 
cality does not need to be invoked (although it may). From this perspective, the violation 
of the Bell inequalities would not constitute a marker of nonlocality. It still remains to be 
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investigated what type of collective or individual phenomena are compatible with the type 
of model introduced in this work. 
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